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We introduce a new fermionic variational wavefunction, suitable for interacting multi-species 
systems and sustaining superfluidity. In this frame, we also introduce a novel way for the treatment 
of entanglement via a new quantum index, which is related to the internal symmetry of the quantum 
state. This index is "ghostly", i.e. it appears only within the quantum state, but it does not appear 
in any observable quantity. Applications include quark matter, nuclei, neutron stars and the high 
temperature superconductors. Spin up and down fermions are, in principle, inequivalent. A wider 
class of Hamiltonians than sheer Bardeen-Cooper-Schrieffer (BCS) type, comprising interaction and 
hybridization between different fermion species, can be treated "exactly", as in the well known 
manner of BCS theory. We present the finite temperature version of the theory, and we discuss the 
appearance of charge and spin density wave order. 



The Bardeen-Cooper-Schrieffer (BCS) wavefunction [l[ |^bcs) = Ylk( u k+ V k c\ ^_ k j_)|0) has set a paradigm for the 
description of fermionic superfiuids. The creation/annihilation operators ct _/cfc )(T describe fermions with momentum 
k and spin a, and |0) is the vacuum state. 

We are interested in a wavefunction for fermionic systems with two or more different species of spinful fermions, 
which can perform a task similar to what I^bcs) can do. 

Let the usual fermionic operators be c\./c x with x = {i, k, cr}, where i denotes the fermion species/flavor. 

We introduce a new quantum index, which is related to the internal symmetry of the quantum state. It serves to 
enumerate the " entangled" components of the quantum state as a function of both momentum and spin. Thereby we 
introduce the new fermionic operators c|. „/c X)/J obeying the anticommutators ({a, b} = ab + ba) 



{ c x,^ii C-y,v\ — , {c x ^, c y,u\ — ^xy $ fiv ■ (1) 

Then, we write the usual c\jc x as the superposition 



= ~7W 51 ' Cx = ~7W ^ Cx ' 5 

v ° 5=1 



(2) 



The usual anticommutation relations of c\/c x are preserved, while 



(3) 



We also introduce 



AU,<5 — U i,k + V i,k c [fc,t,<5 C i,-fc,4.,<5 + S «AT C lfc,t,<5 C j,-fc,i,<5 + s i,k,± c \-k,\.,8 C j,k^,5 • (^) 

A] k £ is a bosonic operator, creating spin singlet pairs of fermions (for triplet pairs c.f. eq. (|60[) below), and 
(i/j)={(M),(2,l)}. 
We form the following multiplet of A[ k s 's 

M l = ^\.k.s=i A,-fc,5=2 A\ k S=z A 2 k S=4 ■ (5) 



This multiplet creates all states with momenta ±fc, so we take N a — 4, as N a such states exist (cf. also after eq. (|56j) 
below). 

The new index allows for the bookkeeping of a superposition of states of a given particle, i.e. same x = {i, k,o~}, 
without the difficulties due to entanglement within the multiplet M k , if the index were removed. In that case, the 
treatment of the coherence factors uj~,v k , Sk is prohibitively complicated. 
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This index is "ghostly", i.e. it appears only within the state |^) below. It does not appear in any 
observable quantity! (C.f. the Faddeev-Popov ghosts in gauge field theories, where the ghosts disappear from 
the final physical outcome.) We note that there is no change whatsoever implied in the Hamiltonian or in the 
representation of any observable. 

Now we introduce the "entangled" state 

m = U M l\o) , (6) 

where the prime implies that k runs over half the momentum space. Note that all A ik s 's in |^) commute with each 
other. 

generalizes I^bcs) and sustains superfluidity. This wavefunction makes sense for two or more fermion species, 
with an interaction between different species. It can obviously be generalized for three or more fermion species. 
Moreover, a similar wavefunction using the new quantum index can be written in the real space representation 
instead of the momentum space one. At the moment, it is not possible to say how close |\&) is to experimental reality. 
It does represent a very promising avenue though, as can be seen from the discussion which follows. 

|\E') is, in principle, relevant for the description of systems such as quark matter [2|, nuclei [3j, |4j, superconducting 
grains [U, cold atoms @, Qi graphene @|, high-temperature superconductors, i.e. both copper oxides @, [loj and iron 
pnictides [HI ], and neutron stars It allows for inequivalence between spin up and down fermions. Plus, it allows for 
the "exact" variational treatment of a wider class of Hamiltonians than sheer BCS type, e.g. comprising interaction 
and hybridization between different fermion species, in the well known manner of the BCS-Gorkov theory [H, 12|. 

The normalization (^l 1 ?) = 1 and the choice of real Ui t k, Vi,ki s i,fc,t> Si,k,l imply 

u i,k + v i,k + s i,k,t + s i,k,l = 1 ; < uf k ,vf kl s1 k ^, sf k ^ < 1 , (7) 

thus allowing to treat these coherence factors as u^k = cos^fc) cos^j^), v^k — sm(8i t k) cos(0j,fc), Si,/c.T = 
cos(<5 ijfc )sin(0 iife ), Si,*,,; = sin((5 ljfe ) sin(0 ijjt ). 

Two fermion species at zero temperature. Calculations are straightforward. E.g. for 2 fermion species with 
dispersions ei.k.o- = U.-k.a we have 

s/N~ ci >fc> t Ml |0) = M c\ _ k i S=1 + s 1At c\ _ k i S=1 ) A\ _ k S=2 A\ k S=3 Al _ k S=i |0) 

-S2,k,i c 2,-fc,l,5=3 A\,k,8=l A\ _ k S=2 ^4-2, -fe, 5=4 |0) • (8) 



Then 



Also, 



(0| M k c\ M c lM Ml |0) = ^{vl k + 4 M + 4 M ) . (9) 



V^Vo c 2M Ml |0) = (v 2tk 4 ^ S=3 + S2M c\ _ k ^ s=3 ) A\ AS=1 A\ _ k S=2 Al _ k S=i |0) 

-«M4 4,-k,l,6=l A \,-k,S=2 4,fc,«=3 A l.-k.S=4 10) , (10) 

which yields 

(0| M k 4 1 fc,t c i.k.t M l 1°) = ~Tr( Vl > k s hk,l + v 2,k s 2 ,k,i) • (11) 

Moreover, 

N C 2 ,_ fc , 4 Ci, fclt Ml |0) = S 1M 4,-fc,5=2 4,fc,«=3 4,-k, ff =4 1°) ~ S 2,k,l 4,M=1 4,-k,«=2 4,- M =4 ' ^ 



and 



(0| Af fe 02,-/5,4.01^^ MjJ |0) = -^-{ui. k si,A,t ~ M 2,fc S2,A,4.) ■ (I 3 ) 
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Using the commutativity of A\ k 5 's and generalizing the previous equations, we obtain ((B) — {9\B\^)) 

n iX a = (cJ tfc , a Ci )fcl0 .> = — {vf >k + si^ a + sl k> _ a ) , (ij) = (1,2), (2, 1) , (14) 
t sgn(cr) 

Zk,a = \ c lk,a C j,k,cr} = (v\M S\M-a + V 2 ,k S 2 ,k,~a) , (15) 

^ o 

gk,a = (c2.-k.-aCl,k,a) = "Tp Sl,fe, CT - U 2 ,k S2,k,-<?) ■ ( 16 ) 

The general Hamiltonian for two fermion species interacting via intra-species potentials V1.2 and via an inter-species 
potential F q , and hybridizing via h^, is 



H = 53 ^> fc > CT C lk,a C ^a + 2 fefc ( c l,fc,<r C 2,fc,o- + 4,fc,<T C l,fe,<x) (17) 
i,k,(T fc,cr 

9 53 ^>9 C !,fc+g,cr C !,p-g,(T' c ^P:0-'Ci,fc : CT + 53 ^9 C l,fe+g,cr C 2 : p-g,cr' c 2,P,o-' Cl,fe, CT , 



i,k,p,q,<T,<r' k,p,q,a,cr' 

with i = 1,2, ^ijfejff = &i,k,(T ~ an d Mi,<7 the chemical potential. Note that the usual BCS pairing potential is just 

the sub-term J2i,k, P v i,k-p c l,fc,t c l,-fc,i c ».-P.i c *.p.t of tne sm gle species potential. 
Considering above, we have ioT{H) = (^\H\^), a Q = 1/N a 

(H) = 53 &,fe,<7 rl ^fc,cr + 2^3 ^fe z fe,cr - - 53 Vi,k-P n i,k,a n^p.cr + 53 Vl -1=° n *-k,<J n i,p,a (18) 
i,k,a k,a i,fc,p,<7 i,k,p,cr 

k~p "Ui.kVi.k u i,pVi,p ^ ' Fk—p z k,rr z p,a + F q= Q S ' 
z,/c,p k,p,tr k, P ,a k.p.a 

with = {(1, 2), (2, 1)}. The first term in the second line is exactly the usual BCS pairing term, and the last term 
is the equivalent inter-species pairing term due to F q . Then, considering the up and down spins as equivalent, we 
minimize (H) with respect to the angles 6^/- and respectively, obtaining 

o = 2c/ l , fe 5^e i ,fc, CT - 2 5I '../.- • ' 51 { v ^=o- v ^-p} n ^+ 2A lk( v ik- u lk) (19) 

<x a p,cr 

-{-^Fq—Q Ui^k / ^ Tlj,p,cr ~H 2 ^ ^ Ffe—p ^,/c,cr ^p,£T 2& V^ft ^ ^ -^fc— p ^i,k,<T ^i,p,cr ^J,p ^j,P,— f) ) 
p,cr p,cr p,£T 

and 

= 2^,^6^-2 J] /ife>W+^X] {^=o-^fe-p} n w + A i,fc sin(2^, fe ) sin(2^ sfe ) (20) 

(T (T p.fj 

— p Si,k,cr (^z,p Si.p.cr ^j-P ^j.p. — cr^} • 

sm( 0j.fc j 

p.(T p,CT p,£7 

with 

p 

= sin(2^, fc ) _ 2g . n2 = ^ cos^) 

and (i, j) = {(1, 2), (2, 1)}. A° fc is the usual BCS gap. Equations (TK)1) - (|2H can be solved numerically (this will appear 

elsewhere). For the number operator n i;PjCT = c\ p a c^ p ^, i d t n itPt „ = [n iiP ^,H] = h p (c\ p ^Cj, p ^ - ct )P)<T Cj )P)<7 ), and 
(d t n i!P>lT ) = 0. 

Our approach is much more general and complete than the works of Moskalenko [13| and of Suhl, Matthias, and 
Walker [14j . which can treat strictly BCS- type inter-species and intra-species interactions only. 
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The finite temperature dependence of the theory can be derived through the equations of motion formalism 
for the Green's functions |l2|, |l5j (d T c x (T) — [H,c x (t)}). We consider the Green's functions 

G l .AP,r-r') = -(T oCi , p , ct (t)4 p ^(t')} , (22) 

G ij>a (p,T-T') = -{T Ci, P)CT (r)ct p)a (r')) , (23) 

fUp,t-t') = (T 4 iCT (r) C t_ p _ CT (r')) , (24) 

4>,t-t') = (T 4 iCT (r) C t_ p _ ff (r')) , (25) 
with = {(1, 2), (2, 1)} and T a denoting imaginary time ordering. We obtain the exact coupled equations 

6(T-T , ) = -(d T + Z i , p , a )G i , tr (p,T-T , )+ ^( T o4 fc _ gff ,(r)c lACT Kr) Cl , p _^(r)4 pff (rO) 

-h p G ji>(T (p,T - t') + F q (T cl k _ qt A^^'(r)ci, P - g Ar)cl P; Ar')) , (26) 

+M^>> T ~ r>) + Y F q {T ct p _, jlT (T)c; )fc+gtff ,(T)c J - lfc , <7 /(r)ci_ 1 , j _ (T (r')) , (27) 

= -(^r + Zj,p,a) G ^a(p,T - T 1 ) + ^ V jtq {To Cj, k - q ,v> ( T ) C ],k,a> (t)c,- (r)cf P)<T (t')) 

-W(p,t-t') + ^ ^<T c] )fc _ ?iCT; (T)c a ,^(T)c i>p _ g;ff (T)cl >P;(T (T')) , (28) 
= (-d T +^ Pt<7 )F} jia (P,r- t') + ]T V 3 , q (To sV g ^(r)4 w (r)c 3 , fc , CT ,(r)4^_ ff (r')) 

+h p Fl a (p,r-T')+ P^cl^ qt Ar)cl^ i Ar)cixAr)4,-p,-a^)) • (29) 

Within the Hartree-Fock-Bogoliubov approximation ((c x c y c z c w ) — (c x c y )(c z c w ) — (c x c z ) (c y c w ) + (c x c w )(c y c z )), and 
after Fourier transforming with respect to r, we obtain the following system of equations at finite temperature T 

1 = (ie„ G x , ff (A;,e„) + Ai Gi, ff (fc, e n ) + B x F^fAr, e„) + C x G 21 , CT (>, e„) + A F^^fc, e„) , (30) 

= {ie n + Ci,fe, CT ) Fl„(k, e n ) + A 2 Gi _ CT (fe, e„) + F 2 i^fc, e„) + G 2 G 2 i,_ CT (fc, e„) + £> 2 ^(fc, e„) , (31) 

- (ie„ - 6,fe,a) G 2 i, ff (fc, e„) + A 3 G lj<7 (fe, e„) + 5 3 Fl _ a (k, e„) + G 3 G 2 i, CT (fc, e„) + £> 3 ^^(fc, e„) , (32) 

= (ie n + F^fe, e„) + A 4 G 1;(T (fe, e n ) + F 4 e„) + G 4 G 21 ,_„(*, e„) + £> 4 F^^fc, e„) , (33) 

and also the equivalent set with the indices 1 and 2 interchanged. The Matsubara energies are e„ = (2n + l)7rT, and 
we suppressed the labels {k, a) from A m , B m , C m ,D m , m = 1 — 4. Here 

M = Y V ^q{ ~ S qA c ip,* Cl >P^ + ^.^'^.P^U-g.a !^-?^)} " F 9=0 Y ( C tp,<?> C 2,P^') > ( 34 ) 
p,q,a' p.cr' 

B\ = - 7J Vj, g (Cl.-fc+g.-CTCl.fe-g.a) , Gl = "/ifc + ^ Fg ( C l , - 9 , <7 4 , fc - <?, CT ) ' ( 35 ) 

-Dl = / ] F q {c2-k+ q -aC\M- q ,a) , 

q 

A 2 = B\ , B 2 = -A, , C 2 = -F q ^{cl k J 2 _ k _ a ) , D 2 = h k + Y p q (<%,k+ q ,A,k +q ,a) > ( 36 ) 

q 

C 3 = Y V2 'l{ ~ S Q>0(4,p,cr> C 2,P,<r') + ^,<T'Sk,p( C 2 ,k- q ,<T C 2,k-q,a)} ~ F q =0 Y ( C hp,<?' Cl 'P'°' ) ' ( 37 ) 
p. q .a' p.cr' 

M = -D 2 , B 3 = Y,F<] (ci,-k+ q ,-<7C2,k-q,cf) ■ D 3 = Y, V2, q ( C 2,-k+q-cfC2,k-q,cr) , (38) 

q q 

A 4 = -F ?=0 <4 , cj^.J ,B A = C X ,C 4 = Dt , D 4 = C 3 . (39) 
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Also, for Ci,k,-a = t,i,k,a, B 3 = -Di and A 4 = -C 2 . 
We have 

M = ^ Vi,g{ - 5g,a ni.p.g + 5 a ^'5 k ,p n 2 ,k- g ,a} - F q=o y^^n 2 ,p,a>, (40) 

p,q,(j' p,a' 

Bi = -a o y^Vi,q m,fc- g ^i,fc- g = A^ fc , (41) 

Cl = -/ife +a ^Fg {Wl,fc-g Sl^-q-a +V 2 ,k-q [Sl,k-q,a + S2,fc-g,-cr] } , (42) 

£>! =^F g5fe _ gi(7 , C 2 = ^,=0 Sk,„ , D 2 = 2h k + C 1 , (43) 

C 3 = V2,g{^g,0 "2,p,o- ~ 5g,<j><ik,p rH,fc-g,<r} + ^=0 ^ "l,p,<7'i (44) 

p,q,cr' p,a' 

D 3 =0, o ^2 V 1,q u 2,k-q v 2 , k - q = -A% k . (45) 
9 

We solve these equations for the case of equivalent up and down spin 5i,fc,-o- = £i,k,a — 6i,fc f° r both fermion species. 
n n ■ n Zi(A;,ie„) ^t/; • \ ^i(fc,ie„) . Z 2 i(k,ie n ) t . Xi 2 (fc,ie„) 

and likewise for G 2 (fc, ie„), i<^(fc, ie„), Gi 2 (fc, ie„) and i*2i(^> * e «)- Here 

= -ie* - 4 (6,fc - ^i) + {Dl - Cl - dD 2 - G 2 B X - 6,fe(2C 3 + &,*)} + ^i[C 3 2 - D§ + &,fc(2C 3 + &,*)] 

-C 2 [d£» 3 + D^Cg + 6, fe )] + 6,fc[CiD 2 - + 2C 3 )] + - Cf) + £>2(£>iA» + CiC 3 ) , (47) 

X! = -Bj e 2 + 2ze n C 2 ^ 2 + D 3 (C 2 - D 2 ) + B^D 2 - Cf - &, fe (2C 3 + , (48) 



Z 21 = -ie n {BiDt + £> 2 (£i,fc + 6,fc + C 3 - Ai) + C 2 D 3 } - £> 3 (BiCi - AiC 2 ) 

-(C 3 + 6,fe)(Si£>i - AiDa) - ei,fe[C 2 ^ 3 + £ 2 (C 3 + &,*)] , (49) 



*i2 = e 2 n C 3 + ie n {A X C 2 - B X C X + C 2 (C 3 + - &,*)} + A»(BiI>i - ^1^2) 

+(<?3 + 6,fc)(SiCi - A X C 2 ) + C 2 (C 2 £>i - CiD 2 + Ci,fe 6,fe) + ^2^3 • (50) 

The denominator is 

^(fc,e) = e 4 + Qe 2 + ,Se + r , (51) 

with Q = B 2 -A 2 -2C 1 D 2 -2C 2 D 1 -C 2 +D 2 HiA2A 1 -£ 1 . k )-&, k (2C 3 H2.k), 5 = 2B 1 (C 2 D 2 -C 1 D 1 )+2D 3 (D 1 D 2 - 
CiC 2 ), and F = - A1) 2 - B 2 ][(6,fc + C 3 ) 2 - £ 2 ] + (dD 2 - C 2 D x f + 2(C 1 £> 2 - C 2 D X ){A X - &, fc )(6,fc + C 3 ). 

Then D(k,E k ) = yields four solutions for the quasiparticle energies E k . Setting K = 27 'S 2 — 72YQ + 2Q 3 , L = 
12Y + Q 2 , N = K - VK 2 - 4L 3 and W = {-2Q + L(2/N) 1 ^ 3 + (N/2) 1 / 3 }/'3 we have 

E(i,2),fc = iVw 7 ±^-VK-2Q-2,S/v / W 7 , £? (3l 4),fc = ± 1 ^-VK - 2Q + 2,5/^ . (52) 

For V U V 2 — > 

£i,fc = -£ 4;fc + 0(yi,T"2) 2 - ^ = -^ + 0(^1,^2) • (53) 
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They depend implicitly on the temperature T through the factors ttj^T), Vi k{T), 8ik,<r{T). The latter can be 
calculated by noting that 

G i>CT (fc,r = 0) = -{cixA,k,a) = T ^G 4 , CT (fc, e „) , (54) 

e n 

FliK t = 0) = (cl ifciCT c] ) _ fei _ CT ) = T £ i^(fc, e „) , (55) 

and through the use of eqs. (|46p. For T < T c , the critical temperature, Fj a (k) become non-zero. 

The superconducting gap in a physical system is probed through various experimental techniques, and it is 
usually extracted from a fitting procedure to some specific theoretical models, including purely phenomenological ones. 
E.g. for the high temperature superconductors these techniques include angle-resolved photoemmission (ARPES), 
tunneling, NMR, Raman scattering etc. As far as the BCS theory is concerned, things are pretty straightforward. In 
our theory, one needs to calculate the precise spectral response in terms of the microscopic parameters of 'F for any 
" gap" -probing experiment, and fit appropriately the data (l6| . 

Charge and spin density wave (CDW/SDW) order can appear in a natural manner, via a simple extension 
of 4 1 . Namely, by allowing the total momentum of pairs to be finite, which is expected to be favored by the finite 
interspecies potential F q . 

For example, we may consider 

A \,k,S = u i,k + V ilk 4,k,\.5 c\_ kxs + Ui,fc,p(fc),t 4,k,t,6 4,P(k),l,S + Vi,- k ,-P(kU 4-k,l,5 4,-P(k),tJ ( 56 ) 
+Si,fc,f 4,k,f,5 c J,-*4,« + s «,*,4. 4,-k,l,6 c ],fe.t,« + **.fc.f > (*),t C lfe,t,<5 c },P{k),l,S + S hk,P(k),l 4-k,l,S c },-P(k),t,S ' 

just for a finite set of momenta k — ±fci, ±&2, ±fcjv- ki in this set is correlated with some other kj — P{k i ) 1 
and fcj = P(kj), (i,j=l-N), and so are their opposite momenta — fcj and — kj. To fix the ideas further, considering 
a 2-D system with tetragonal symmetry C4, we may take N=4, and the correlated pairs as (fci, £2), (^3) ^4) & n d 
(— k\, — fe), ( — A?3, — fei). In this case A Q = 8(= 4 + 4). In principle, the new index could be a continuous variable, 
instead of an integer, if a continuous range of momenta would be correlated with a given k. This seems to be the 
actual physical case. We might consider in addition the correlations (k\, — k 2 ), (— k\, fc 2 ), (k 3 , — £4), (~fc 3 , £4) etc. 
For a pair of such correlated momenta (k,p) we obtain 

(Cj ^ Ci,p,<r) — 7^ { Sgn((j) Si,k.p.~a ^z,p Si,p,k, — a ~t~ ^i.k.p. — a ~l~ ^i,p,k,— <j] 

~t~^J)fe,p, — cr ^J,fe, — cr "i" Sj,p,k, — <T ^j^p, — cr } 7 (^^) 

( c i,-fc,(T c j,p,<r) = T^ - { — ^ s S n ( cr ) [^i.fc Si t k,p,~a + fj,p Sj'.p.fc.-er] + s j,k,-a Vj t k,p, — a + s i,p, — <J v i,p,k,-a} 1 (58) 

(4,-fe.a C i,P,-<r) = (4,-fe.a C j\P,-^ = • ( 59 ) 

We may further introduce spin-triplet pairing terms such as 

^i.fe.CT c i,fe, CT ,5 C i,P(fc),CT,<5 + *i,fe,o- C !,fc,o%<5 C j,P(/c),er,<5 ' (^) 

etc. in At . 5 , which also contribute to CDW/SDW order, yielding 

(4-k,a C i,P,-a) = Jj-{ S 3,k,-a *J,*,p,-ff + s j,p,cr *j,p,k,<7 + Sgn(cr) Wt.p.k.c - V ijk Wi,k,p,-v] } , (61) 



1 

A„ 



{sj,fe,-(r Wj,k,p,-<r + s, lP , CT Wi iPt k,o + sgn(cr)[uj )P i,\ p ,fc, CT - tt,*, Pl - CT ]} ■ (62) 



We note that m principle it is possible to have non-zero ex pec tation values for charge and spin density 
( c l fc+Q o- c i,fe,±o-) - with both i = j and i ^ j - for some Q-range 17|, while the anomalous propagators = 0. 
This regime could correspond to the pseudogap phase of the copper oxide superconductors flij |. 

In OUil (treating different models though!) the coexistence of charge and spin density wave order with supercon- 
ductivity was explored. We emphasize that, as far as we can currently see, this coexistence is not compulsory, though 
possible, in our approach. 



The relevance of the BCS states to the calculation of T c for the multilayer copper oxide superconductors [10( 
provided a motivation for this work. 

The author is indebted to Gregory Psaltakis and John Smyrnakis for valuable discussions. Alex Gezerlis pointed 
torefs. gi. 

* E-mail address : kast@iesl.forth.gr , giwkast@gmail.com 
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